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Abstract 

Multielectron bubbles (MEBs) differ from gas-filled bubbles in that it is the Coulomb repulsion 
of a nanometer thin layer of electrons that forces the bubble open rather than the pressure of 
an enclosed gas. We analyze the implosion of MEBs subjected to a pressure step, and find that 
despite the difference in the underlying processes the collapse dynamics is similar to that of gas-filled 
bubbles. When the MEB collapses, the electrons inside it undergo strong accelerations, leading to 
the emission of radiation. This type of sonoluminescence does not involve heating and ionisation 
of any gas inside the bubble. We investigate the conditions necessary to obtain sonoluminescence 
from multielectron bubbles and calculate the power spectrum of the emitted radiation. 
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I. INTRODUCTION 



Multielectron bubbles (MEBs) are typically micron-sized cavities in liquid helium, con- 
taining a nanometer thin layer of electrons on the inside of the surface of the bubble. MEBs 
were first observed in experiments studying electrons on the surface of liquid helium beneath 
which an anode was present^. When the charge density of the 2D-layer of electrons exceeds 
a critical value, a surface instability sets in^ and a bubble forms, carrying up to 10 8 electrons 
to the anode. The electron system inside such a MEB is of particular interest, since it forms 
a spherical two-dimensional electron gas whose density can be tuned over a wid e range of 
values depending on the number of electrons in the bubble and the pressure applied on the 
helium^. 

Unlike in gas-filled bubbles, the force preventing the bubble from imploding is not due 
to the pressure of a gas inside the bubble, but due to the Coulomb repulsion between the 
electrons in the MEB. At zero external pressure, the radius of the bubble is determined 
by the balance between the Coulomb repulsion and the surface tension of the heliuirA When 
the pressure p on the helium is nonzero, the bubble radius satisfies the equation 

e 2 N 2 

where e is the electron charge, N is the number of electrons in the MEB, e is the dielectric 
constant of helium, and a is the surface tension of helium. In the regime pR^jo ^> 1 the 
bubble radius scales as p -1 / 4 . In this paper we investigate what happens when a bubble 
with radius Rb(p = 0) in equilibrium at zero pressure is suddenly subjected to an increase in 
pressure. Its radius is larger than the equilibrium radius at pressure p, R\,{p = 0) > R\,{p) 
and we expect the bubble to shrink. In section II we derive and investigate the equations of 
motion for the bubble dynamics, and in section III we apply these equations to analyze the 
collapse dynamics of the multielectron bubble. 

The fact that inside MEB's there are electrons (charged particles) instead of neutral gas 
atoms necessarily has another interesting consequence: the electrons in the compressing 
bubble are subject to accelerations, and therefore they generate electromagnetic radiation. 
The breathing mode of the bubble gives an estimate of the response time of the helium 
surface - which typically is microseconds for a 10000 electron MEB. In that time span, the 
bubble radius shrinks a distance of typically a micrometer. The velocity corresponding to 
this is of the order of micro meter /microsecond or m/s. But the acceleration this generates is 
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on the order of a micrometer /microsecond 2 which is quite large (10 6 m/s 2 ). As the number 
of electrons increases, the initial bubble radius increases and the maximum acceleration 
increases as well. We will show that ultrashort light pulses can be produced by a sudden 
pressure increase in the helium surrounding the MEB. 

Explanations of sonoluminescence in gas-filled bubbles^^ have relied strongly on the 
emissivity of the heated gas in the collapsing gas-filled bubble^, as opposed to the mechanism 
of sonoluminescence proposed here for multielectron bubbles. Nevertheless gas-filled bubbles 
are known to get to temperatures where at least weak ionisation of the gas must occur—, 
and electrons can emit light, being accelerated through scattering from the neutral atoms. 
Although emission from accelerated electrons can contribute to the sonoluminescence of gas- 
filled bubbles^, it's precise role is still disputed^. For MEBs we propose that accelerating 
electrons form the principal mechanism of sonoluminescence. Moreover the light emission is 
expected to be coherent and coming from the bubble surface, whereas in gas-filled bubbles 
the emission is incoherent black-body emission from the interior of the bubble. 



II. EQUATIONS OF MOTION FOR THE BUBBLE SURFACE 

The bubble surface at time t can be described by the function R(6, 0, t) that gives the 
distance from the center of the bubble to the surface, in the direction specified by the 
spherical coordinates 8, <fi. This function can be developed in spherical harmonics 

oo £ 

R(0, 0, t) = R h (t) + a ^(t)Ye,m(^ 0), (2) 

i=\ m=-l 

where i?b is the angle-averaged radius of the bubble. The coefficients ae m describe the 
deformation of the bubble from its spherical shape. To derive the equations of motion for 
the bubble surface, we start from the bubble Lagrangian derived iiA 

ae,m\ 2 , (3) 

where p represents the mass density of liquid helium. For small amplitude deformations of a 
spherical bubble, we know that a^ m {t) = a^ m (0) x exp{iuiet} where the ue are the spherical 
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ripplon frequencies 3 . These frequencies can be obtained from (jSJ) by putting Rh(t) = R eq 
where R eq is the equilibrium radius: 



7 + 1 



eq 



o{t 2 + i + 2) + 2pR, 



eq 



e 2 iV 2 



(4) 



In the Lagrangian we have neglected the influence of the redistribution of the electrons 
along a deformed bubble surface. This effect can be described by introducing the (spherical) 
plasma oscillation modes of the electron system, and calculating the coupling between plas- 
mons and ripplon^. In a recent study of this coupling Klimin et a/.— found that although 
a bubble deformation does lead to a redistribution of charge along the surface, the resulting 
shift of the ripplon frequencies is small. The reason for the weakness of the ripplon-plasmon 
coupling is the difference in frequency of these modes: the ripplon frequencies typically lie 
in the MHz-GHz range, whereas the plasmon frequencies lie in the THz range. 
The equation of motion for the radial component is 

d dC dC 



dt dR h dR h 



(5) 



which yields 
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The first line corresponds to the Rayleigh equation (with an additional electronic term) 
describing the collapse of spherical bubbles^. The second line in © contains additional 
terms due to the deformation of the bubble. 

The equation for the deformation components {ae }in (t)} is obtained from 

d dC dC 



dt da 



da 



Lm 



which leads to 



di,m + o — ai )Tn — — j- 



a(f + £ + 2) + 2pR h 



e 2 N 2 ( 



(7) 



(8) 



This equation is similar to the Plesset-Prosperetti equation^ in the theory of gas-filled 
bubbles in fluids. The set of equations © and (JBJ need to be solved given appropriate 
initial conditions. 
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Salomaa and Williams 14 studied small-amplitude oscillations of multielectron bubbles 
using a similar set of equations. They neglected the contribution of the deformation am- 
plitudes {a^ m }'s in their equation for the radial component Rb(t). This approximation 
is valid if ag jTn <C i?b- They found that when the (time averaged) bubble radius satisfies 
R cq < [N 2 e 2 /(IQ-K(re)} 1 ^, the t = 2 deformation can grow indefinitely, suggesting an insta- 
bility of the bubble. In subsequent work^ it was found that the I = 2 ripplon frequency can 
indeed vanish, corroborating the results of Salomaa and Williams. However, an extension of 
ref.— for large amplitude oscillation (i.e. beyond the regime a^ m <C R^) showed that there 
exists a metastability barrier preventing the bubbles from fissioning^. 



III. COLLAPSE OF A MULTIELECTRON BUBBLE 



Having obtained the equations of motion for the bubble shape, we can proceed to solve 
them for the particular case of an MEB subjected to a sudden increase in pressure. We start 
with a bubble whose angle-averaged radius i?b is equal to the equilibrium radius R eq (po) at 
a given pressure p of the liquid helium surrounding the MEB. This bubble can be either 
spherical (W, m : a^ m (t = 0) = 0) or deformed (3£, m : a^ m {t = 0) ^ 0). At time t = the 
pressure of the liquid helium jumps from p Q to a higher value p. The angle-averaged radius 
is then larger than the equilibrium radius at pressure p and the bubble will collapse. 



A. Spherical bubble 



First we consider the collapse of an undeformed bubble. Then all a^ m 's are zero and the 
dynamics are governed by the Rayleigh equation 1 - 1 

2 
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2 \R^ K = ~p~R 
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2aR h + pR\ 



(9) 



in which an additional term due to the Coulomb repulsion of the electrons appears. This 
equation can be integrated once by multiplying left and right hand side by 2R^R h , resulting 
in 

' e 2 N 2 



Jt ( RlRl 



1 d 
2iip dt 
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\2eR h 



+ AnaRl + —pR 3 h 



(10) 
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so that 
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[U(R h ) - U(R b (0))} 



(11) 
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where U(R) 
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This can again be integrated to yield 
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[U(R h (0))-U(r)} 



dr. 



(13) 



Rb(0) 



Note that the solution R(t) is periodic, with the period given by T = t(R tnrn ) where -Rtum 
is the classical turning point in the potential U(r). In figure ^ the time evolution of the 
radius of a spherical N = 10000 electron bubble subjected to a sudden increase in pressure 
is shown. The top panel shows the pressure increase, and the subsequent panels show the 
radius, the radial velocity and the radial acceleration as a function of time divided by the 
period T. The initial conditions are -Rb(O) = and -Rb(O) = 1.06441 /im - this corresponds 
to the N = 10000 electron bubble in equilibrium at zero pressure. The periods are T = 0.795 
/xs, 0.570 /is and 0.462 fjs for pressure steps of 1,2, and 3 kPa, respectively. 

The radial acceleration shows a pronounced peak at the time when the MEB radius 
is minimal (the 'time of collapse'). Upon increasing the pressure, the maximum radial 
acceleration grows and the oscillation of the radius becomes more strongly anharmonic. 
There is no dissipation term in our equation for the radius, so that the radius bounces back 
to its original value during the periodic oscillation. Dissipation will reduce this value and 
will let the bubble perform a damped oscillation around its new equilibrium radius R eq (p). 
As can be seen from figure 1, the time evolution of the radius of the MEB is very similar 
to that of a gas-filled bubble subjected to a shock wave^. In the latter case, the restoring 
force opposing the collapse is the increase in pressure of the gas trapped in the bubble. In 
the present case, the bubble contains no gas (except helium atoms at vapor pressure which 
is negligible at low temperatures), and it is the increased Coulomb repulsion that gives rise 
to the restoring force. 

B. Deformed bubble 

For small-amplitude deformations, a£. m <C R\>, we can neglect the \a^ m \ terms in the 
equation (0) for the bubble radius, and use the solution for R^it) discussed in the previous 
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section directly in the equation (JHJ) for the deformation amplitudes. This equation can then 
be written as 

a t , m (t) + \{t)a itm (t) + u 2 {t)a^ m {t) = 0, (14) 

with the 'drag coefficient': 

\{t) = 3R h {t)/R h {t), (15) 
and the time dependent ripplon frequency 

£+1 



P Rl(t) 



e 2 N 2 

a(£ 2 + £ + 2) + 2pR h (t) - 



(16) 



47reRl{t) 

The time dependence of these coefficients (|T5 |l . (fT5|l appears through the time-dependence 
in i?b- The differential equation (JTljl can be solved numerically. The main question of this 
subsection concerns whether the deformations will vanish, grow without bound, or stay finite 
during the collapse of the bubble. From the numerical analysis we find that there are stable 
and unstable modes of deformation. 

Note that even at a fixed pressure there are stable and unstable modes, as mentioned 
at the end of section II. Moreover, at higher pressure, more modes are unstable. When a 
pressure step Ap is applied, a mode that was stable at the initial pressure p may be unstable 
at the final pressure po + Ap. 

An example of the time evolution of the amplitude of a stable mode is shown in figure 
El For this figure, we start with an N = 10 4 electron bubble at po = that has an 
£ = 25 ripplon present. In this example, £ is chosen arbitrarily, but sufficiently high so that 
we(po+Ap) > 0. Under the equilibrium conditions, the deformation amplitude oscillates with 
the corresponding ripplon frequency. When the pressure is suddenly increased by Ap = 10 3 
Pa, we find that the oscillation perseveres, but that the oscillation amplitude is modulated 
with the same periodicity as the bubble radius. When the bubble radius is minimal, the 
oscillation amplitude of the deformation mode is maximal. This can be understood by 
substituting the ansatz a^(t) = A(t)e luJt in the equation of motion (JHJ). Collecting the 
imaginary parts of this equation, we arrive at 

Ait) R h , , 

from which 

A{t) = CR- 3/2 (t). (18) 
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The constant of integration is given by C = A(0)R h (0). This solution is shown in figure El 
by the dashed curve. It fits well with the envelope of the full numerical solution (full curve), 
even though we started with a rather large deformation amplitude (10% of the radius) in 
order to get a clear figure. For stable modes, we thus find that during the collapse of a 
deformed bubble the deformations do not vanish. The deformation amplitude grows, but is 
bounded at all times. 

An example of the time evolution of an unstable mode is shown in figure 01 As in figure 
121 the initial state is that of an MEB with N = 10 4 electrons at pressure zero, subjected 
to a pressure increase of 10 3 Pa. The different curves in figure 01 represent different initial 
conditions, as listed in the legend. We chose the £ = 2 mode because it is known to have a 
vanishing frequency at positive pressures and thus may lead to instabilities. The numerical 
solution of (fPf|) shows that the deformation amplitude a 2 does indeed grow to a point where 
it becomes comparable to the bubble radius. This indicates that the collapsing MEB in 
helium becomes unstable and may break up, in contrast to gas-filled bubbles in cryogenic 
liquids 1 -. Nevertheless, at the time of smallest radius, the deformation amplitude is still 
smaller than the bubble radius. The point of instability (a 2 = -Rb) is reached after the 
collapse. So, although the bubble may not survive long after the first collapse, it will reach 
the regime where the radial acceleration is maximal one time. How long the bubble survives 
after it has reached its smallest radius, depends strongly on the initial conditions of the 
deformation. To study the exact dynamics of the instability (the possible fissioning), the 
full set of equations ©,© needs to be solved. This lies beyond the scope of this paper as 
we want to investigate the collapse dynamics and luminescence during collapse, but not the 
possible breaking up of bubbles during the explosion that follows after the implosion. 



IV. LIGHT RADIATION FROM AN IMPLODING MEB 



A. Incoherent emission from a spherical bubble 



The electric radiation field by an accelerated charge is given by 



E(r,R,t) = 4 nX[nX ^ )] 



c 2 |r-R(t)| 



(19) 



ret 
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where R(i) gives the trajectory of the charge q, r is the point in which we want to evaluate 
the electric radiation field, n is the unit vector along r — R(i), and c is the speed of light in 
vacuum. Note that there is a non-radiative part of the electric field (of order [r — R(t)]~ 2 and 
higher). This is usually neglected at large distances from the accelerated electron. However, 
if the radiative field is cancelled, for example by symmetry considerations, then these higher 
order fields come into play so that the Maxwell equations remain satisfied. The subscript 
"ret" means that the electric field at time t is generated by the accelerating electron at a 
time t' earlier by the delay it takes the light to reach r. In our situation, this retardation 
is of no consequence. From figure El it can be seen that though the radial velocity of the 
electrons is fairly high, it is not relativistic (the acceleration, be it large, only takes place in 
a very short time scale). The power radiated per unit of solid angle per electron is 

with 9 the angle between the direction of the acceleration and n. Thus the maximum 
radiated power is 

We first consider a light pulse that might be emitted by a single electron in an MEB. The light 
pulse intensity is peaked at the instant when the bubble reaches its smallest radius, stops 



R 



is largest. 



imploding and starts expanding again. At this point, the radial acceleration 
This is similar to what happens in sonoluminescence in gas-filled bubbles. For bubbles with 
10000 electrons, and a pressure step of 3 kPa, a pulse of roughly 50 ns is emitted, as can be 
seen from figure ^ bearing in mind that in that case T = 0.462 /is. Increasing the pressure 
step will shorten the length of the pulse, and increase the maximum acceleration reached 
(see figure Q). In figure 0] we show the width At pu i se at half maximum of the peak in the 
radial acceleration, as a function of pressure for different numbers of electrons. Large MEBs 
(N = 10 7 ), such as those created in the experiments reported in Refs.— , have a very high 
compressibility resulting in ultrashort pulses and very large radial accelerations. Figure El 
shows the power as a function of time (and frequency, in the inset) for an extreme case: an 
N = 10 7 electron bubble subjected to a pressure step of 101.3 kPa. The ultrashort pulse 
reaches a maximum power of 100 juW (i.e. 10 -11 W per electron). Of course, one can ask 
whether for such an extreme case, the equation of motion @ remains valid up to the time 
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that the smallest radius is reached. If the surrounding helium is heated up by the radiation, 
or if the electrons penetrate into the helium, then the bubble implosion may have different 
dynamics. Nevertheless, it seems safe to expect that equation (jHJ) remains valid for a large 
portion of the collapse, and that for not so extreme cases the qualitative results hold. These 
results are that if the number of electrons or the pressure step is increased, the pulse duration 
shortens, the emitted spectrum shifts to higher frequencies, and the maximum power of the 
pulse increases. 



B. Sonoluminescence from deformed bubbles 



The results of the previous subsection were for one accelerating electron, or for a collection 
of electrons emitting radiation incoherently. If the electric field of all the electrons is added 
coherently, then for an undeformed bubble the radiation of the collapsing shell of electrons 
will vanish due to the spherical symmetry, as we shall prove. 

We investigate the case where one of the deformation amplitudes a^ m is nonzero. In 
studying modes of vibration where m/0, the restriction a^_ m = a| m applies because the 
surface R(fl, t) has to be described by a real function, and thus 

R(fi,t) = {R h (t) + 2Re[a e>m (t)Y £ , m (fl)]}n n (22) 
= R b {t) [1 + a e , m (t)P™ (cos 6) cos(m0)] n n , (23) 

where P™ is the associated Legendre function, n$] is the unit vector normal to the surface 
and the rescaled deformation amplitude ag m is given by 



We evaluate the electric field at a point 

r = r (e x sin a sin (3 + e y sin a cos (3 + e z cos a) , (25) 

where we assume that r ^> R, i.e. the radiation field is evaluated at a large distance from 
the bubble. Furthermore, we will assume that de )7n -C -Rb, which is valid when the bubble is 
near its minimum radius, and when the initial deformation satisfies a^ m (t = 0) <C Rb(t = 0). 
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Substitution of in (JTHJ) results in 



^ NeR h (t) 
E(r > i?) « — ^ x 



iVe R h (t) 



R h (t) v 7 i2 b (t) 
dn [1 + a Am (t)P;"(cos0) cos(m0)] [0(0, 0)e r - e R ] 



c 2 rR h (t) 

where the angular function is given by 

(0, <ft) = cos a cos + sin a cos j3 sin 6* cos + sin a sin /3 sin sin < 
The components of the electric field are, for r 3> -R 



(26) 
(27) 



(28) 
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\ E >) 



q a e>m (t)R h (t) 
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cos mc 
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^ 0(0, 0) sin a cos (3 — sin cos 
(0, 0) sin a sin /3 — sin sin 
0(0, d>) cos a — cos0 



(29) 



From the integration over the <fi coordinate, it is clear that only the m = 0, ±1 terms result 
in a nonzero radiation field. The integration over the angle finally yields 



E y 

\ E -) 



^Ne R h (t) 
T^~rR h (t) 



t 
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cos a sin a cos /3)ai i0 



sin 2 a cos 2 /3)ai t 



\ 



(cos a sin a sin (3) ai >0 — (sin 2 a sin /3 cos 
(sin 2 a)ai j0 — (sin a cos a cos (3)ai t i 



(30) 



Note that only the i = 1 mode contributes to the radiation field. This proves our earlier 
statement that the radiation fields of electrons on an accelerating spherical shell cancel out. 
The emitted power is then given by 



P(a,P) = N 2 P 



' sin a cos (3 cos a ^ 



sin a sin (3 cos a 



LV 



sin 2 a 



Qi,o 



' sin 2 a cos 2 /? — 1 
sin 2 a cos /3 sin (3 
cos a sin a cos /3 



(31) 
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with P = (4n/3)e 2 Rl/c 3 . Using typical values iV = 10 6 ,a^ m ph 0.01 Rb, po = and 
Ap = 100 kPa, we find a maximum power of the order of 100 fiW (for a pulse length of the 
order of 10 fs). Although this may seem a lot, because of the short pulse time only a small 
amount of energy is emitted, corresponding to 10-100 photons with frequency 1/(10 fs). 
This will make the detection of the light difficult in practice. Photons of successive collapse- 
expansion cycles may be collected, but dissipation may quickly dampen the oscillations. 
If the bubble oscillations after a single pressure step dampen out quickly, repeating the 
pressure step using an ultrasonic wave of e.g. 1 KHz can increase the number of photons to 
be collected. Another possibility may be to work with larger bubbles since the total number 
of electrons affects very strongly the maximum peak power, through the coherency factor 
N 2 in (|31|) and the strong dependence of on N. For example, a N = 10 5 electron bubble 
subjected to the same Ap = 100 kPa pressure step has a maximum peak power of 10 -2 
Z^VV (for ci£ m ~ 0.1 i?b). In the extreme case of large bubbles (N > 10 7 ) subjected to large 
pressure steps (Ap of the order of 10 5 Pa), the power radiated can become large enough 
to affect the collapse dynamics significantly, and expression (pTT|) is no longer valid. Two 
factors will act to reduce the maximum power: incoherency effects (such as those caused by 
the finite time needed for a pressure wave to travel across the bubble) and the smallness of 
the deformations (typically, we assumed ai m -C .Rb)- 

The last factor in (J3"T| is the geometrical factor giving the angle-dependence of the ra- 
diation and relating its strength to the deformation amplitude. The collapsing deformed 
bubble will radiate anisotropically. Figure El shows the angular dependence of the radiated 
power. The (£,m) = (1,0) mode will emit most strongly in the equatorial plane of the 
bubble, whereas the (£,m) = (1, ±1) modes will emit most strongly along meridians. 

V. DISCUSSION AND CONCLUSION 

When an MEB is subjected to a sudden increase in pressure p$ — > po + Ap, its radius is 
larger than the new equilibrium radius R eq (po) > R e q(Po + Ap) and it shrinks to a final radius 
that is smaller than the new equilibrium radius R eq (po + Ap), after which it expands again. 
This oscillatory motion of collapses and expansions occurs on a time scale of microseconds, 
but the oscillation becomes greatly anharmonic, in that most of the reduction in radius 
occurs in a small fraction of the total cycle. This leads to large accelerations when the 
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radius of the bubble is smallest. When deformations are present, the normal modes of 
deformation with a ripplon frequency larger than the inverse time scale of the collapse are 
stable. Modes with small ripplon frequencies (such as the £ = 2 quadrupole oscillation 
modes) grow unbounded and can lead to instabilities and fissioning of the collapsing bubble. 
However, we find that at the time when the bubble has reached its smallest radius, and the 
radial acceleration is maximal, the deformation mode is not yet diverging. This conclusion 
is important for the possibility to observe sonoluminescence of the collapsing MEBs. In a 
purely spherical bubble, coherent emission of radiation by the accelerated electrons cancel 
out due to symmetry requirements. Thus, to produce a radiation field, the collapsing bubble 
should be deformed. We find that in a collapsing bubble where the i = 1 mode of deformation 
is present, the radiation fields of the electrons do not cancel. One might wonder how this 
particular mode can be excited so that sonoluminescence can be generated in an experiment. 
It turns out that it is not necessary to artificially excite this mode: Plesseti^ showed that gas- 
filled bubbles that implode near a surface naturally develop a deformation corresponding 
exactly to the I = 1 mode, and we expect that this surface proximity effect will also be 
present for imploding multielectron bubbles. 

In the derivation, it was assumed that the temperature is low enough so that the vapor 
pressure of helium inside the MEB is negligible. It is worth noting that if helium vapor is 
present inside the bubble, the pressure of this vapor will increase as the bubble collapses 
and a mist of helium may be formed that can affect the dynamics of the bubble collapse^. 

In summary, we have investigated the dynamics of collapsing multielectron bubbles in 
liquid helium and the possibility that imploding MEBs emit radiation. We find that these 
bubbles behave in a similar manner as gas-filled bubbles, even though the fundamental 
mechanism involved in the bubble dynamics is different. In MEBs the force counteracting 
the collapse bubble is Coulomb repulsion, whereas in gas-filled bubbles this force is due to 
the gas pressure inside the bubble. The collapse of a MEB subjected to a pressure step leads 
to the emission of radiation due to the acceleration of the shell of electrons at the bubble 
surface, whereas sonoluminescence in gas-filled bubbles is due to a more complex combined 
effect of heated gas emissivity and bremsstrahlung. 
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FIG. 1: A multielectron bubble with N = 10 4 electrons at zero external pressure, is subjected at 
time t = to a pressure increase of Ap = 1, 2, 3 kPa (full, dashed and dotted curves, respectively). 
In this figure, the panels from top to bottom show the pressure, the bubble radius, the radial 
velocity and the radial acceleration as a function of time. All these quantities are cyclic with 
period T = 0.795 ps, 0.570 fj,s and 0.462 /xs for pressure steps of 1,2, and 3 kPa, respectively. 
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FIG. 2: The time evolution of the amplitude of deformation associated with the I = 25 ripplon 
is shown (full curve) for aniV = 10 4 electron bubble subjected to a pressure increase of 1 kPa at 
time t = 0. The dashed curve shows the analytical estimate for the the envelope of the oscillating 
deformation amplitude. The dotted curve depicts the radius of the bubble. 
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FIG. 3: The time evolution of an unstable (£ = 2) mode of an N = 10 4 bubble subjected to a 
pressure increase of p = 10 3 kPa. The deformation amplitude grows larger than the bubble radius 
(so that the bubble is unstable) some time after collapse (i.e. after the bubble radius reaches its 
minimum), depending on the initial conditions. 
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FIG. 4: The width at half maximum A£ pu i sc of the peak in the radial acceleration of the collapsing 
bubble as it reaches its smallest radius is shown as a function of the pressure step causing the 
collapse, for different values of the number of electrons in the MEB. In the inset, the maximum 
power of the emitted radiation pulse is shown as a function of the pressure step and the number 
of electrons. 
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Time (fs) 

FIG. 5: For very large MEBs (10 7 electrons) subjected to a large pressure step (100 kPa), the 
emitted pulse of radiation calculated with (R)]).(|19)) can become extremely short and reach 100 fiW 
peak power. This also results in a broad spectrum of the emitted radiation, as shown in the upper 
inset. The lower inset shows the bubble radius as a function of time (cf fig. 1). 
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FIG. 6: The maximum power (expression (|31|0 of the radiation from a collapsing deformed bubble 
is shown as a function of the spherical angles and <j). The power is expressed relative to maximum 
power as a function of the angles. The top two figures are for the a^o ^ 0, ai 5 i = case and the 
bottom two figures for the ai.o = 0, ai i 7^ case. The illustrations on the right show the regions 
on the sphere (with the z-axis in the vertical direction) where radiation is emitted (shaded regions) 
in both cases. 
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